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ABSTRACT 

We present an analysis of the nnclear infrared (IR, 1.6 to 18 /im) emission of the 
ultraluminous IR galaxy UGC 5101 to derive the properties of its active galactic nu¬ 
cleus (AGN) and its obscuring material. We use new mid-IR high angular resolution 
(0.3 — 0.5 arcsec) imaging using the Si-2 filter (Ac = 8.7 ^m) and 7.5 — 13 /tm spec¬ 
troscopy taken with CanariCam (CC) on the 10.4m Gran Telescopio CANARIAS. We 
also use archival HST/NICMOS and Subaru/COMICS imaging and S'piteer/IRS spec¬ 
troscopy. We estimate the near- and mid-IR unresolved nuclear emission by modelling 
the imaging data with GALFIT. We decompose the Spitzer/IKS and CC spectra 
using a power-law component, which represents the emission due to dust heated by 
the AGN, and a starburst component, both affected by foreground extinction. We 
model the resulting unresolved near- and mid-IR, and the starburst subtracted CC 
spectrum with the CLUMPY torus models of Nenkova et al. The derived geometrical 
properties of the torus, including the large covering factor and the high foreground 
extinction needed to reproduce the deep 9.7 /rm silicate feature, are consistent with 
the lack of strong AGN signatures in the optical. We derive an AGN bolometric lumi¬ 
nosity Lboi ^ 1.9 X 10"*® erg s“^ that is in good agreement with other estimates in the 
literature. 
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1 INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs) are among the 
most luminous objects in the local universe, with both 
their luminosities (Lir > lO^^L© emerging mainly in the 
far-infrared) and space densities being similar to those of 
quasars (e.g. Sanders & Mirabel 1996). ULIRGs have large 
molecular gas concentrations in their central kpc regions 
(e.g. Downes & Solomon 1998) with densities compara¬ 
ble to those of the stars of elliptical galaxies. For increas¬ 
ing infrared (IR) luminosities, their mid-IR slopes and sil¬ 
icate depths also increase, suggesting that their nuclei are 
more obscured and more compact (Stierwalt et al. 2013 and 
references therein). Kormendy & Sanders (1992) proposed 
that ULIRGs evolve into elliptical galaxies through merger- 
induced dissipative collapse. In this scenario, these mergers 
first go through a luminous starburst phase, followed by a 
dust-enshrouded active galactic nucleus (AGN) phase, and 
finally evolve into an optically bright, naked quasar once 
they either consume or shed their shells of gas and dust 
(Sanders et al. 1988a). 

The target of interest in this paper is the nearby (z = 
0.039; Kim 1995) ULIRG UGC 5101 (IRAS 09320-f6134) 
with a luminosity distance d = 168 Mpc (1 arcsec = 755 
pc) for Ho = 73 km s“^Mpc“^, fim = 0.27, I2a = 0.73, 
and an IR luminosity Lir = 10^^ L© (Sanders et al. 1988a). 
In the optical and near-IR the galaxy has a single, very 
red nucleus within a disturbed morphology, suggestive of a 
recent interaction and merger (Sanders et al. 1988b; Surace 
et al. 2000; Scoville et al. 2000). 

Optically the AGN has been classified as a Seyfert 1.5 
(Sanders et al. 1988a), a LINER (Veilleux et al. 1995), 
and a Seyfert 2 (Yuan et al. 2010). Mid-IR spectroscopy 
taken with ISO (Genzel et al. 1998) as well as ground-based 
L-band spectroscopy (Imanishi et al. 2001) demonstrated 
the presence of a powerful nuclear starburst. Using ground- 
based high angular resolution mid-IR imaging, Soifer et al. 
(2000) showed that « 60% of the IRAS flux at 12 /rm comes 
from the central 4 arcsec (~ 3 kpc) and that almost half of 
this arises in an unresolved core. Based on its IRAS col¬ 
ors, UGG 5101 is classified as a cold, starburst-dominated, 
far-IR source. However, Chandra and XMM-Newton data 
detected the Fe-Ka line, which together with the low ob¬ 
served I/2-iokev/LiR and I/2-iokev/L[oiiil ratios likely in¬ 
dicate the presence of a Compton-thick AGN (Ptak et al. 
2003; Imanishi et al. 2003; Gonzalez-Martin et al. 2009). 
At radio wavelengths, the nucleus of UGG 5101 has a high 
brightness temperature (T > 10^ K) with a morphology re¬ 
sembling that of a jet, again pointing to the presence of an 
AGN (Lonsdale et al. 1995, 2003). Using Spitzer IR spec¬ 
trograph (IRS, Houck et al. 2004) data, Armus et al. (2004, 
2007) detected the mid-IR fine structure line of [Nev] at 
14.3 pm, which confirms the presence of an AGN in the nu¬ 
cleus of this galaxy. 

UGC 5101 offers an ideal scenario for studying the in¬ 
terplay between nuclear activity and star formation in a 
nearby ULIRG. The ^5—10 arcsec spatial resolution of 
Spitzer/IRS only allowed to study this interplay in ULIRGs 
on kpc-scales, although with extremely good sensitivity. For 
instance, Veilleux et al. (2009) used different indicators on 
the Spitzer/IRS spectrum of this galaxy and estimated an 
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Figure 1. The left panel is the fully reduced GTC/CC 8.7 /rm (Si-2 filter) image of UGC 5101, in the middle is the PSF standard star 
shown in a linear scale. The right panel is the PSF-scaled (at the 65% peak level, see section 3.1.2) subtracted image of UGC 5101. 
The lowest contour is 3(t over the background, with the next contours in 2a steps (except in the PSF image). The orientation of all the 
images is North up, East to the left. The physical scale of the UGC 5101 images is 755 pc/". 


average IR luminosity, whereas Armus et al. (2007) esti¬ 
mated a 30% contribution from a starburst component. 

Here we present new high spatial resolution (~ 
O.Sarcsec ~ 230 pc) mid-IR imaging and spectroscopy of 
UGC 5101 obtained with CanariCam (CC, Telesco et al. 
2003; Packham et al. 2005) on the 10.4 m Gran Telescopic 
CANARIAS (GTC). We use these observations to disentan¬ 
gle the unresolved mid-IR emission, which is presumably 
due to dust heated by the AGN, from that due to star for¬ 
mation in the nuclear region of the host galaxy. By mod¬ 
elling the unresolved near-IR and mid-IR emission with the 
CLUMPY torus models of Nenkova et al. (2008a, b) we aim 
to derive the physical and geometrical properties of the ob¬ 
scuring material surrounding the AGN of UGC 5101. The 
paper is organized as follows. Section 2 presents the obser¬ 
vations and data reduction. In Section 3 we estimate the 
unresolved emission in the near-IR and mid-IR using the 
imaging and spectroscopic data. In Section 4 we model the 
unresolved near and mid-IR emission of UGC 5101 and de¬ 
rive the properties of the AGN and its obscuring material. 
Section 5 gives the conclusions of this work. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 GTC/CC mid-IR imaging and spectroscopy 

As part of an ESO/GTC large programme (ID: 182.B-2005, 
PI: A. Alonso-Herrero) to observe a large sample of local 
AGN, in January 6th, 2014, we observed in queue mode 
UGC 5101 with CC on the GTC in El Roque de los Mucha- 
chos Observatory (La Palma, Spain). CC uses a Raytheon 
320 X 240 Si:As detector that covers a field of view (FOV) 
of 26 arcsec x 19 arcsec on the sky with a plate scale of 
0.0798arcsec pixeU^. We obtained an imaging observation 
with the Si-2 filter (Ac = 8.7 /rm and width AAcut = 1.1 
^m at 50% cut-on/off) using the standard mid-IR chop-nod 


technique. The chop and nod throws were 15 arcsec. The 
airmass during the observation was 1.2. We observed the 
mid-IR spectro-photometric standard star HD 79354 (Cohen 
1999) before the science target. Allowing for overheads, this 
was 20 min before the science observations. The standard 
star is used to perform the photometric calibration, esti¬ 
mate the full width half maximum (FWHM) of the science 
observation, and perform the point spread function (PSF) 
subtraction. We measured the FWHM of the data by fitting 
the radial profile of the standard star with a Gaussian func¬ 
tion and obtained a value of 0.4 arcsec. Table [T] summarizes 
the details of the observations. 

We also obtained long-slit mid-IR spectroscopy of the 
nuclear region of UGC 5101 on the same night. We used 
the low spectral resolution 10 ^m grating, which covers the 
A-band ~ 7.5 — 13 ^m with a nominal spectral resolution 
of R = A/AA ~ 175, and a 0.52 arcsec width slit. In Table[T] 
we list the position angles (PA) of the slit, measured from 
the North to the East. The observing sequence was as fol¬ 
lows. We first took an acquisition image of the galaxy with 
the Si-2 filter, then centered the slit on the peak mid-IR 
brightness, and finally integrated for the on-source integra¬ 
tion times given in Table [T] We used the same chop-nod pa¬ 
rameters as for the imaging observations. We also observed 
standard stars using the same observing sequence to provide 
the photometric calibration, the telluric correction, and the 
slit loss correction. We also used the acquisition image of 
the standard star to obtain the FWHM of the spectroscopic 
observation (see Table [T]). 

We reduced the CC data using the RedCan pipeline for 
the reduction and analysis of ground-based mid-IR CC and 
T-ReCS imaging and spectroscopic data (Gonzalez-Martin 
et al. 2013). Briefly, the reduction process of the imaging 
data includes sky subtraction, stacking of the individual im¬ 
ages, and rejection of bad images. The flux calibration of the 
galaxy image is done using the observation of the standard 
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Table 1. Summary of the GTC/CC observations. 


UGC 5101 

Imaging 

Spectroscopy 

Date 

2014.01.06 

2014.01.06 

NOB x ton (s) 

1 x 1224 

1 x 1242 

Standard 

HD 79354 

HD 79354 

NOB X Gn-STD (s) 

1 x 60 

1 X 60 

FWHM (arcsec) 

0.4 

0.3 

PA (degree) 

0 

90 


Notes.— Nos is the number of repetitions performed, ton is 
the on-source integration time of UGC 5101 and ton— STD is the 
on-source integration time of the standard star. The FWHM is 
measured from the image of the standard star. The PA value for 
imaging is the orientation of the detector on the sky, while for 
spectroscopy it indicates the orientation of the slit. 


star. Fig. [T] (left panel) shows the fully-reduced 8.7 fj,m CC 
image of UGC 5101 together with that of the standard star 
(middle panel). As also found by Soifer et al. (2000) at 12 /tm 
using the Keck telescope, UGC 5101 is clearly extended in 
the mid-IR (see Section 3.2). 

For the spectroscopy the first three steps of the data re¬ 
duction are the same as for the imaging. Additionally, Red- 
Can performs the two-dimensional wavelength calibration of 
the target and standard star spectra using sky lines. Then, 
the trace determination is made using the observation of the 
standard star. The last steps of the data reduction include 
the extraction of the spectra either as point sources or ex¬ 
tended sources and finally, the correction for slit losses in 
the case of point source extractions. 


2.2 Archival HST/NICMOS and 
Subaru/COMICS imaging 

Our target was observed in a single orbit in November 
7 of 1997 with the Hubble Space Telescope (HST) using 
camera 2 of NICMOS (Thompson et al. 1998) as part of 
the Guaranteed Time Observations. This camera uses a 
256 X 256 HgCdTe array with pixel scales of 0.0762 arcsec 
and 0.0755 arcsec pixeP^ in x and y, respectively, provid¬ 
ing a ~ 19.5 arcsec x 19.3 arcsec FOV. We use the im¬ 
ages of UGC 5101 taken with the broad-band F160W fil¬ 
ter (Ac = 1.60/rm) and the medium-band F222M filter 
(Ac = 2.22 ^m). We downloaded the five reduced and cal¬ 
ibrated images in each of the F160W and F222M filters 
from the Mikulski Archive for Space Telescopes (MAST). 
We combined the individual exposures of each filter using 
the average of the images. These images were originally pre¬ 
sented and analysed in Scoville et al. (2000). 

We also downloaded a fully reduced and calibrated im¬ 
age of UGG 5101 taken in the Q-band (Ac = 17.7 /rm) from 
the German Virtual Observatory (GAVO) hosted at the 
Virtual Observatory (VO). The image was obtained with 
COMICS on the Subaru Telescope in Manna Kea, Hawaii. 
It has a FOV and a pixel scale of 42 arcsec x 32 arcsec and 
0.13 arcsec pixel”respectively. For full details on the Q- 
band observations and data reduction we refer the reader to 
Asmus et al. (2014). 


2.3 Spitzer/IRS spectroscopy 

UGC 5101 was observed with Spitzer/IRS using the short- 
low (SL) and long-low (LL) modules covering the follow¬ 
ing spectral ranges: SLl; 7.4 — 14.5/rm, SL2; 5.2 — 7.7 pm, 
LLl; and 19.5 — 38.0 pm, LL2; 14.0 — 21.3/rm, with a spec¬ 
tral resolution of 7? ~ 60 — 120. The SL slit width is 
3.7 arcsec, whereas that of the LL module is 10.5 arcsec. 
We downloaded the fully-calibrated staring mode spectrum 
from the Cornell Atlas of Spitzer/IRS Source (CASSIS v6; 
Lebouteiller et al. 2011). CASSIS has identified this object 
as point-like and therefore it uses its optimal extraction rou¬ 
tine to produce the best flux-calibrated spectrum and ensure 
the best S/N ratio. We stitched the different module spectra 
together by using the shortest module SL2 flux as reference 
spectrum for the other modules. The scaling factors were 1.1 
for LLl and LL2. 


3 ANALYSIS 

3.1 Unresolved nuclear emission from the imaging 
data 

Even at high spatial resolutions, in the near-IR up to A ~ 2 
pm, extended stellar emission arising in the host galaxy con¬ 
taminates and even dominates the nuclear fluxes of AGN, es¬ 
pecially in type 2s. At longer wavelengths A > 3 pm contam¬ 
ination by stellar photospheric emission is greatly reduced 
(Alonso-Herrero et al. 2001 and references therein). Our goal 
in this section is to use the high angular resolution imaging 
data to isolate as much as possible the unresolved nuclear 
emission before we attempt to infer the properties of the 
obscuring material and the AGN luminosity of UGG 5101 
(Section 4). To do so, we need to remove any extended nu¬ 
clear emission not directly related to the dusty torus and/or 
emission arising from dust heated by young massive stars 
(see Alonso-Herrero et al. 2011 and references therein). 

Although in the past there have been some attempts 
to measure the unresolved nuclear emission of UGC 5101 at 
near- and mid-IR wavelengths (e.g. Soifer et al. 2000; Scov¬ 
ille et al. 2000; Haan et al. 2011; Imanishi et al. 2014; Asmus 
et al. 2014), the reported values may still be contaminated 
by stellar emission and emission due to nuclear star forma¬ 
tion. Therefore, we measure the unresolved emission consis¬ 
tently in the archival HST/NICMOS H- and A-band images 
and COMICS Q band images as well as in the new CC im¬ 
age. To do so, we use the two-dimensional fitting algorithm 
GALFIT (Peng et al. 2002) to derive both the unresolved 
nuclear emission and extended emission, as GALFIT allows 
to simultaneously fit an arbitrary number of components 
(e.g. PSF, Sersic profile, exponential disk, and/or Gaussian 
functions) to a galaxy image. It is possible, however, that 
even at these high angular resolutions (0.15 — 0.40”) there 
is still some contamination from other components (stellar 
emission and star formation, see Section 4). 

3.1.1 Near-IR emission 

To model the unresolved nuclear emission we generated the¬ 
oretical Tiny Tim PSFs for the NIGMOS camera 2 F160W 
and F222M filters. The Tiny Tim PSFs have only three free 
parameters, namely the central coordinates {xc and y^) and 
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Figure 2. Upper panels: to the left is the combined HST/NICMOS 2.22 fim image (F222M filter) of the nuclear region of UGC 5101, 
in the middle the GALFIT model image (see Table 2 for the parameters), and to the right the residual image after subtracting the 
GALFIT model from the original image. The images are shown in a linear scale, rotated so that North is up. East is to the left. The 
lowest contour is 3cr above the background, with the next contours in 2a steps. Lower panels: same as upper panels but for the reduced 
GTC/CC 8.7 iim (Si-2 filter) image. The physical scale of the the images is 755 pc/”. 


the total flux. The total unresolved emission is computed 
by integrating the flux over the PSF image assuming that 
it contains 100% of the light. To fit the extended surface 
brightness distribution, that is, the host galaxy emission, we 
used a Sersic profile (see Peng et al. 2002 for more details 
on profiles). 

To constrain the unresolved emission in the near-IR and 
its uncertainty we use different combinations of profiles that 
always include a PSF profile to represent the unresolved 
emission. These combinations of profiles are: a PSF and a 
Sersic profile with all the parameters free, a PSF plus a 
Sersic profile with fixed n = 1, a PSF plus a Sersic pro¬ 
file with fixed n = 4, and finally one PSF plus two Sersic 
profiles including a Gaussian profil^il with all the parame- 

^ A Gaussian profile is a special case of the S&sic profile when 
n = 0.5 and the size parameter is the FWHM instead of r^g. 


ters left to vary freely. We find that the models that include 
a PSF plus two Sersic profiles with free parameters and a 
PSF plus one Sersic profile with n=4 produce fits with sim¬ 
ilar ~ 0.35 and that the unresolved emission measured 
on the NICMOS images did not change by more than 1% 
at both wavelengths. We then conclude that the unresolved 
component is well constrained at both near-IR wavelengths. 
Fig-[2]shows the combined NICMOS 2.22 fim image together 
with the GALFIT model image and the residual image after 
subtracting the GALFIT model from the image. 

Since the fits with different models are of similar quality, 
we estimate the uncertainty of the parameters in Table 2 as 
the standard deviation of the values given by all models. 
These uncertainties are larger than those estimated for a 
single model. 

The unresolved emission detected in UGC 5101 con¬ 
tributes 2% and 17% of the continuum emission at 1.60 
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Table 2. Results from the GALFIT modelling using a PSF + Sersic profile model. 


Instrument 

A 

(pm) 

./unresol 

(mjy) 

/Sersic 

(mJy) 

n 

reft 

(pc) 

a/h 

PA 

(degree) 

NICMOS 

1.60 

< 1.3 

21 ±6 

4* 

1753 ± 728 

0.60 ±0.03 

82 ±2 

NICMOS 

2.22 

5 ± 2 

25 ±6 

4* 

1091 ± 584 

0.60 ±0.03 

82 ±3 

CC 

8.70 

36 ±4 

169 ± 10 

2.0 ±0.2 

1127 ±65 

0.60 ±0.02 

78 ± 1 


Notes. — Column 1: instrument. Column 2: central wavelength of the filter. Column 3: unresolved nuclear emission. Column 4: 
emission of the Sersic component. Column 5: index of the Sersic profile. *For the 1.60 and 2.22 pm model the index was fixed. Column 
6: effective radius in pc of the Sersic component. Column 7: axis ratio of the Sersic component. Column 8: PA of the major axis of the 

S&sic component measured East of North. 


and 2.22 pm, respectively, measured over the /fST/NICMOS 
Camera 2 FOV. Our unresolved fluxes at H and K band are 
lower than those reported by Scoville et al. (2000) using 
the same NICMOS/HST images. This is because they only 
subtracted the underlying galaxy as estimated from the ad¬ 
jacent pixels instead of modelling the host galaxy emission. 
The H band image was also analyzed by Haan et al. (2011) 
using a similar decomposition technique as used here, where 
they found similar results for PSF contribution and Sersic 
parameters, within the uncertainties. We adopt our own val¬ 
ues for the analysis that follows, since our estimation of PSF 
uncertainties follows the same methodology as that used to 
measure the unresolved component at K and mid-IR bands, 
and it considers a larger range of models than in the analy¬ 
sis by Haan et al. The host galaxy parameters are similar to 
those found in previous studies at both near-IR bands (e.g. 
Scoville et al. 2000, Haan et al. 2011, Medling et al. 2014). 

3.1.2 Mid-IR emission 

As done for the HST/NICMOS images, we also run GAL¬ 
FIT on the CC 8.7 pm image to estimate the unresolved 
and integrated emission of the galaxy. For the unresolved 
component we used the CC observation of the standard star 
taken just before the galaxy observation. We assume that the 
FWHM did not vary on those time scales and therefore the 
emission of the standard star can be used to represent the 
unresolved emission of UGC 5101. We fit a PSF and Sersic 
component with all parameters left free to vary. The best 
fit had a reduced ~ 1- The uncertainties for the param¬ 
eters are the standard deviation of the parameters derived 
by considering different convolution boxes. The lower pan¬ 
els of Fig. [2] show the original CC 8.7 pm image together 
with the 2D GALFIT model and the residual image. Af¬ 
ter subtracting the GALFIT model we obtained a residual 
of less than 1.5% when compared to the flux in the origi¬ 
nal image. Table |3] lists the unresolved and Sersic fluxes at 
8.7 pm. The unresolved component contributes ~ 12% of the 
observed emission at this wavelength. The uncertainties in 
the reported profile fluxes are the formal photometric errors 
reported by GALFIT added quadratically to the standard 
deviation (6% for both the Sersic and the PSF components) 
of the fluxes obtained by considering different convolution 
box sizes in GALFIT and to the 10% flux calibration un¬ 
certainty (see Table [si . 

We also estimated the unresolved and integrated emis¬ 
sion at 8.7 pm with two independent methods routinely used 
in the literature for ground-based mid-IR imaging. The first 
one is referred to as PSF-scaling photometry and is useful 


only in sources that are dominated by the nuclear emis¬ 
sion (see e.g. Mason et al. 2012). In this method, we mea¬ 
sured the flux within a circular aperture of 0.5-arcsec radius. 
Next, we scaled the PSF-star image to the centroid of the 
galaxy emission and measured the flux within the same aper¬ 
ture around the scaled PSF star. This method provides the 
maximum contribution from an unresolved nuclear source. 
For UGC 5101 we obtained the integrated emission using a 
large aperture with a 9.2 arcsec radius, which is the radius 
at which the curve of growth flattens. 

The second method we term PSF-subtraction (e.g. 
Radomski et al. 2002, 2003; Soifer et al. 2000; Levenson et 
al. 2009; Ramos Almeida et al. 2009, 2011; Garcfa-Bernete 
et al. 2014). In this method, we first matched the PSF-star 
image (see Fig. [1] right panel) to the peak of the galaxy 
emission, that is, at a 100% level. Then we subtracted the 
scaled standard PSF-star from the galaxy image at differ¬ 
ent percentage peak levels until we obtained a fiat profile 
in the residual image at < 0.3 arcsec. For the CC 8.7 pm 
image this was achieved at a subtraction level of 60% of the 
peak intensity. We then measured the unresolved component 
by integrating the emission in a 9.2 arcsec-radius aperture 
on the scaled PSF-star image. We computed the integrated 
flux by adding the unresolved component to the extended 
component as measured by integrating the flux in the same 
aperture radius on the galaxy image after subtracting the 
scaled PSF image (see Fig. [1] right panel). 

We list in Table [3] the unresolved fluxes at 8.7 pm us¬ 
ing the PSF-scaling and PSF-subtraction methods. We com¬ 
puted the uncertainties adding in quadrature the photomet¬ 
ric erroi0 and the uncertainty in the flux calibration (10%). 
For the PSF-scaling method we assumed a 6% uncertainty 
due to time-variability, estimated from the variation in the 
signal from the standard stars, and a 13% uncertainty due 
to PSF-variations. These uncertainties are derived in Ma¬ 
son et al. (2012) from several standard stars observed dur¬ 
ing the same night at mid-IR with T-ReCS and Michelle 
on the Gemini North telescope. The uncertainties assumed 
were added in quadrature to the errors previously described. 
For the PSF-subtraction method we added in quadrature, in 


2 Calculated as -I- 'rback'^pixZ-^pix-ring, where 

is the number of pixels inside the aperture considered, A^pix—ring 
is the number of pixels inside a ring of 80 pixels width around 
the source used to estimate the background level and its standard 
deviation (Tback (Reach et al. 2005). The background in our images 
is completely fiat, and this makes the second term of the error 
equation almost negligible. 
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Table 3. Unresolved and integrated emission at 8.7 fim using three different methods. 


Method 

,/unresol 

(mJy) 

R apert 

(arcsec) 

fint 

(mJy) 

R apert 

(arcsec) 

GALFIT 

36 ±4 

9.2 

205 ± 11 

NA 

PSF-scaling 

40 ±7 

0.5 

240 ± 41 

9.2 

PSF-subtraction 

41 ± 12 

9.2 

168 ± 32 

9.2 


Notes.— NA=No assigned. 

addition to all the above error contributions, a 7% of uncer¬ 
tainty due to ambiguity in the best subtraction of the un¬ 
resolved source. The unresolved fluxes estimated with the 
PSF methods are in good agreement with the value esti¬ 
mated with GALFIT. 

In order to assess the robustness of decomposition and 
whether there might be a dependency with the adopted PSF, 
we use the acquisition image of the standard star for spec¬ 
troscopy obtained 2 hours after the science image (see Table 
[TJ and repeat the calculations of GALFIT decomposition of 
UGC 5101 with this PSF model. We find an unresolved flux 
27±6 mjy, which is consistent with the value derived before. 
We also find that the Sersic parameters are consistent with 
those previously derived and listed in Table 

We also used GALFIT to measure the unresolved flux 
in the Subaru/COMICS Q-band image at 17.7 fim. In this 
case we used a Gaussian profile with a FWHM ~ 0.6 — 
O.Sarcsec, which is the range FWHM of the standard star 
(see Asmus et al. 2014), to represent the unresolved com¬ 
ponent and a Sersic profile for the extended component. 
We considered three Sersic profiles: one free, one with fixed 
n = 4, and one with fixed n = 1. The best fit was obtained 
with a FWHM = 0.65 arcsec for the unresolved component 
and n = 2.1 ± 1.8. This produced an unresolved flux at 
17.7 fim of 158 ± 32 mJy. This is lower than the unresolved 
flux calculated by Asmus et al. (2014) using the same data 
set but consistent with that in Soifer et al. (2000) using 
Keck imaging. The error in the unresolved emission is the 
quadratic sum of the standard deviation of the flux mea¬ 
sured by fitting different profiles and the formal error from 
GALFIT. 


3.2 Starburst subtracted nuclear emission from 
the spectroscopic data 

Both the Spitzer/IRS and the GTC/CC spectra show emis¬ 
sion from polycyclic aromatic hydrocarbon (PAH) features 
(see Figs 4 and 5, and Armus et al. 2004). This indicates 
the presence of a component due to star formation activ¬ 
ity in the regions of ~ 393 pc and 2.9 kpc probed by 
CC and IRS, respectively. In Section 4 we will also fit the 
starburst-subtracted nuclear GTG/CG spectrum to infer the 
properties of the obscuring material, presumably the dusty 
torus, around the AGN of UGC 5101. Before we do so, we 
need to remove the emission due to star formation. Spectral 
decomposition in the mid-IR has proven to be an efficient 
method to disentangle the AGN emission from that due to 
star formation activity in the host galaxy in luminous in¬ 
frared galaxies (LIRGs), ULIRGs, AGN, and sub-millimeter 
galaxies (e.g. Pope et al. 2008; Valiante et al. 2009; Nardini 
et al. 2010; Mullaney et al. 2011; Alonso-Herrero et al. 2012; 


Table 4. Results from the spectral decomposition of the IRS and 
CC spectra. 


Parameter 

Spitzer/IKS 

GTC/CC 

Power-law index C 2 

2.0 ± 0.2 

2 . 0 (fixed) 

Power-law rxg .7 

5.1 ± 0.6 

6.4 ± 2.0 

Starburst Txg .7 

0.9 ± 0.4 

0.3 ± 0.2 

8.7pm /sBsubs (mJy) 

79± 17 

53 ±35 

12 pm /sBsubs (mJy) 

62± 18 

42± 15 

Reduced 

2.3 

3.0 


Notes.—/sBsubs is the flux obtained from subtracting the ob¬ 
served spectrum and the best-fit starburst model at 8.7 fim and 
12 fim, respectively. 


Ramos Almeida et al. 2014b; Hernan-Caballero et al. 2015, 
Garcfa-Bernete et al. 2015). 

For the spectral decomposition we used the model pro¬ 
posed by Pope et al. (2008). This model includes three main 
components: (1) a starburst dominated by the broad fea¬ 
tures due to PAH molecules up to ~ 18 pm (Puget & Leger 
1989; Allamandola et al. 1999), (2) a power-law (or warm 
blackbody) continuuirQ that represents the obscuring mate¬ 
rial around the AGN (Genzel et al. 1998), and (3) extinction. 
The starburst and power-law components can be affected by 
different degrees of extinction, which is modelled as a fore¬ 
ground dust screen. The model is expressed as: 

F, = -h C4/.(PAH)e""''"", (1) 


where Ti, is the extinction law, ca and cs are the power-law 
component and starburst extinctions, C 2 is the index of the 
power-law, /jy(PAH) the starburst component, and ci and 
C 4 are the normalizing factors. 

We used a total of 24 starburst galaxies from the lit¬ 
erature including IR-bright galaxies, LIRGs, and ULIRGs 
(Sturm et al. 2000; Brandi et al. 2006; Rieke et al. 2009; 
Sargsyan et al. 2011). We used the extinction laws of Ghiar 
& Tielens (2006) for the local interstellar medium (ISM) 
and the Galactic Center (GC), as they include the absorp¬ 
tion due to silicates around 10 and 18 pm. We performed 
a minimization and fitted simultaneously the five free 
parameters ci, C 2 , cs, C 4 , and cs for each of the starburst 
templates and extinction laws. 
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Rest-frame Wavelength (/xm) 


Figure 3. Upper panel: spectral decomposition of the Spitzer/IRS ~ 5 — 35/rm spectrum from the central 4arcsec (~ 2.9 kpc) of 
UGC 5101. The black line is the observed spectrum with the grey shaded region representing the Icr uncertainty. The red line is the 
best-fit model (the pink shaded area is the Itr uncertainty of the fit), which is the sum of the fitted starburst template (blue line) and the 
power-law continuum component (green line). Both components are assumed to be obscured by different levels of foreground extinction 
using the Chiar &: Tielens (2006) law for the local ISM (see text and Equation 2). Bottom panel: residuals of the fit, where percentage 
is defined as observed spectrum minus modelled spectrum, divided by modelled spectrum. 


3.2.1 Spectral decomposition of the Spitzer/IRS spectrum 

Considering that the local ISM extinction law of Chiar 
& Tielens (2006) covers a range from 1.2 to 27 to¬ 
gether with the high S/N and spectral range covered 
by the Spitzer/IRS spectrum 5 — 35 fim compared with 
the GTC/CC spectrum, we started by decomposing the 
Spitzer/IRS spectrum between 5 — 26 /rm. Then, in Sec¬ 
tion 3.2.2, we use the starburst galaxy obtained from the 
best Spitzer/IRS spectral decomposition fit to perform the 
spectral decomposition of the nuclear CC spectrum. The 
best fit for the Spitzer/IRS spectrum is achieved using the 
LIRG ESO 244-G012 spectrum (blue line) from Sargsyan 
et al. (2011) and a power-law component with an index of 
C 2 = 2.0 ± 0.2 (green line). The index of the power-law is 
similar to that found in the mid-IR for Seyfert Is and Seyfert 
2s (Armus et al. 2007; Honig et al. 2010; Ramos Almeida et 
al. 2009, 2011). The fitted power law component in the mid- 
IR is consistent with the presence of a hot dust component 
(T > 300 K) found by Armus et al. (2007) for this galaxy 
from the fitting of the infrared SED of this galaxy. 

We also found that the local ISM Chiar & Tielens (2006) 
extinction law produced a better fit to the data than GC one. 
As can be seen from Fig. 4, the fit reproduces well the mid- 
IR emission except for the bluest part of the spectrum. This 
is due to the presence in UGC 5101 of strong absorption due 
to water ice under the 6.2 pm PAH feature (see Spoon et al 
2002 and Armus et al. 2004), which is not present in the 
extinction law. Adding the water ice absorption using the 
laboratory analogs of Gerakines et al. (1995) only improved 


® Although the use of a power-law continuum in the mid-IR (up 
to A ~ 35 pm) is not physically motivated, it represents well the 
AGN heated dust continuum which usually peaks at A ~ 30 — 
50;rm (see e.g. Alonso-Herrero et al. 2011). 


marginally the fit in this spectral region and the power-law 
component remained the same. 

TableUsummarizes the relevant parameters of the spec¬ 
tral decomposition of the IRS spectrum. The errors in the 
derived parameters take into account the 10% uncertainty 
assumed for the Spitzer/IRS spectrum and the standard de¬ 
viation of the parameters by considering different initial val¬ 
ues for the index of the power-law C 2 (see Section 3.2.2). The 
starburst-subtracted fluxes (observed, not corrected for ex¬ 
tinction) at 8.7 pm and 12 pm (observed wavelengths) con¬ 
tribute proximately 54 and 50% of the observed emission 
in the Spitzer/IRS spectrum, respectively. The starburst- 
subtracted fractional contribution at 12 pm is in agreement 
with that estimated by Hernan-Caballero et al. (2015) us¬ 
ing the same IRS spectrum but an independent method. At 
rest-frame wavelengths longer than approximately 15 pm the 
power-law component appears to dominate the continuum 
emission in the IRS spectrum, although at Arest > 35 pm 
the starburst component has an almost equal contribution. 
This is also in agreement with the 30% starburst contri¬ 
bution found by Armus et al. (2007) for UGC 5101 using 
different mid-IR diagnostics. 

Using the foreground extinction and silicate optical 
depth ratio Av/rxg.r = 18 from Roche & Aitken (1984) for 
the local ISM the fitted foreground extinction affecting the 
power-law component is Av = 90 ± 11 mag, which is neces¬ 
sary to reproduce the deep 9.7 pm silicate feature. This high 
value is also consistent with the presence of dust lanes in 
the nuclear region, as seen in the HST optical image of this 
galaxy (Surace et al. 2000). 

Finally, we note that using more typical starburst 
galaxy spectra, like those of Brandi et al. 2006 we obtained 
similar starburst-subtracted fluxes at 8.7 and 12 ^m, extinc¬ 
tions at 9.7 pm and power-law index to those obtained with 
the galaxy ESO 244-G012, although the x2 values of the 
decomposition are larger. 
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3.2.2 Spectral decomposition of the GTC/CC spectrum 

To determine the extraction type of the GTC/CC spec¬ 
trum (point source versus extended source) we measured 
the nuclear flux on the CC 8.7 pm image within a 1 arcsec- 
diameter circular aperture. We obtained a flux of 57 ± 5 
mjy (after correction for point source emission). A compar¬ 
ison with the unresolved emission at 8.7 pm (Table 3) shows 
that for the CC slit width we can approximate the emis¬ 
sion within the slit as a point source. For the extraction 
as point source RedCan uses an extraction aperture that 
increases with wavelength to account for the decreasing an¬ 
gular resolution. It also performs a correction to account for 
slit losses. We Anally resampled the nuclear CC spectrum to 
60 data points to improve the S/N ratio and for later fit¬ 
ting with CLUMPY torus models (Section 4.2). The errors 
are computed as the quadratic sum of the 10% photometric 
uncertainty and the rms of the four data point rebinning. 

We did the spectral decomposition of the GTC/CC nu¬ 
clear spectrum as for Spitzer/YRS. However, because of the 
limited spectral range of the CC spectrum, we fixed the star- 
burst template to the one that provided the best fit to the 
IRS spectrum (i.e. ESO 244-G012) and used the Chiar & 
Tielens (2006) extinction law of the local ISM. We found 
that the best-fit with a reduced = 2.4 is achieved with 
an index of the power-law C 2 = 4.0. Although this value 
is consistent with the range of spectral indices found by 
Buchanan et al. (2006) in a sample of Seyfert 1, 2, LINERs 
and starburst galaxies, it is larger than that obtained from 
the Spitzer/YRS> decomposition (see Table 14)1. Therefore, we 
study the behavior of the C 2 index in the Spitzer/YRS> and 
CC decomposition by comparing the resulting parameters 
when we vary C 2 between 0 and 4. From this analysis we 
found that the power-law and starburst extinctions and, the 
starburst-subtracted fluxes at 8.7 and 12 pm are well con¬ 
strained (see Table 4). Thus, we chose to fix the power-law 
index to C 2 = 2.0 for the spectral decomposition of the CC 
spectrum for consistency. This value agrees well with that 
derive from the spectral decomposition of the Spitzer/YRS 
spectrum and those measured in Seyfert 1 and 2 nuclei in 
the mid-IR (Honig et al. 2010; Ramos Almeida et al. 2009, 
2011 ). 

These results are not highly dependent on the starburst 
template library adopted. If for instance, we fit the full spec¬ 
tral library of Brand et al. (2006), we find an equally valid 
starburst decomposition (x^ = 3.0) for NGC1222, and sim¬ 
ilar values for extinction, and starburst-subtracted fluxes. 

We show the best-fit model (dashed red line and shaded 
region) to the observed GTC/CC nuclear spectrum (black 
line and grey shaded region) in Fig. |4l with the green and 
blue lines showing the individual contributions from the 
power-law continuum and the starburst components, respec¬ 
tively. As can be seen from this figure, the best-fit model re¬ 
produces well the data and in particular the emission of the 
11.3 pm PAH feature. We therefore subtracted the starburst 
component from the CC spectrum, so in the next section we 
can model the starburst-subtracted mid-IR emission. 

We list in Table |4] the relevant parameters of the fit to 
the CC nuclear spectrum. The errors of the parameters are 
estimated as the standard deviation of the parameters given 
by different models that fix the value of C 2 between 0 and 
4. The fitted foreground extinction affecting the power-law 



Rest-frame Wavelength (^m) 

Figure 4. Upper panel: spectral decomposition of the GTC/CC 
nuclear (inner 393 pc) spectrum of UGC 5101. We resampled the 
observed spectrum to 60 data points (see text). Lines are as in 
Fig- El where the power-law index is fixed to C 2 = 2.0. The filled 
red dot is the 8.7 pm unresolved flux estimated from the GALFIT 
modelling in Section 3.1.2., whereas the orange filled dot is the 
12 pm unresolved flux from Soifer et al. (2000). Bottom panel: 
residuals as calculated in Fig. 4. The hatched region represents 
approximately the spectral range of low atmospheric transmis¬ 
sion. 

continuum from the spectral decomposition of the CC nu¬ 
clear spectrum Ay = 115 ± 36 mag is higher than that from 
the IRS spectrum, although consistent within the uncertain¬ 
ties, possibly indicating that on nuclear scales the heating 
source is more embedded and the continuum source is ex¬ 
tended (see below). 

The starburst-subtracted fluxes, not corrected for ex¬ 
tinction, at 8.7 pm and 12 pm (observed wavelengths) from 
the decomposition of the CC spectrum (see Table 4) con¬ 
tribute approximately 67% and 59% to the nuclear CC spec¬ 
trum, respectively. Within the uncertainties these values 
are in good agreement with those estimated by Hernan- 
Caballero et al. (2015) for the AGN component. The 
starburst-subtracted fluxes at 8.7 pm and 12 pm from the 
Spitzer/YYlSi and CC decomposition are consistent with each 
other, within the uncertainties (see Table 4). In general 
for nearby Seyfert nuclei the spectral decomposition of IRS 
spectra is able to recover both the shape and flux of the AGN 
component measured from high angular resolution ground- 
based mid-IR spectroscopy (see e.g. Ramos Almeida et al. 
2014b; Hernan-Caballero et al. 2014; Garcfa-Bernete et al. 
2014). In the case of UGC 5101, as it is more distant than 
these nearby Seyferts, it is possible that the IRS spectrum 
has a strong contribution from extended dust components 
which are more difficult to disentangle from the dust heated 
by the AGN. 

We And, however that the 8.7 pm starburst-subtracted 
nuclear flux from the decomposition of the GTC/CC nuclear 
spectrum is in excellent agreement to within the uncertain¬ 
ties with the unresolved emission derived from the GALFIT 
model of the CC imaging data (Section 3.1.2 and Table |3]). 
Soifer et al. (2000) estimated that at observed A = 12.5 pm 
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Figure 5. Cartoon showing the geometry of the CLUMPY dusty 
torus models of Nenkova et al. (2008a, b). According to this model 
each dusty cloud has an optical depth ry, and the clouds are 
distributed between the inner radius R^, determined by the dust 
sublimation temperature, and the outer radius Ro = YR^, with 
Y a free parameter. The angular distribution of the clouds is a. 
The torus viewing angle from the observer is i. 

the core component contributed approximately 50% of the 
flux within a 4 arcsec diameter aperture 185 ± lOmJy). 
Within the uncertainties, their inferred unresolved flux is 
consistent with the starburst-subtracted GTC/CC spectrum 
at the same wavelength (~ 91 mjy) but not with the fitted 
power-law continuum (see also Fig. 5). This result shows 
that at this part of the spectrum the power-law is under 
estimated while the starburst component is well fitted. 

4 ESTIMATING THE PROPERTIES OF THE 
AGN AND THE OBSCURING MATERIAL 

4.1 The CLUMPY torus models 

The CLUMPY model (Nenkova et al. 2002, 2008a, b) holds 
that the dust surrounding the central engine of an ACN is 
distributed in clumps. An illustration of this model is shown 
in Fig. [5] The clumps have a radial distribution oc with 
a radial extent Y = Ro/Rd, where Ro and Rd are the outer 
and inner radius of the toroidal distribution, respectively. 
The inner radius is defined by the dust sublimation temper¬ 
ature, which is assumed to be T^oh ~ 1500 K, and the ACN 
luminosity as: Rd — 0.4 (1500 K rg“^)^ ®(LAGN/10'‘®erg 
g-i)0.5 pg model each clump has the same opti¬ 

cal depth {tv). The number of clouds along the line of sight 
(LOS) at a viewing angle i is Alos(*) = 

where (Ttorus is the angular width occupied by the cloud dis¬ 
tribution and No is the average number of clouds along a ra¬ 
dial equatorial ray. The CLUMPY models have been proven 
to reproduce the nuclear IR emission of Seyfert galaxies and 
quasars (e.g. Ramos Almeida et al. 2009, 2011, 2014a,b; 
Nikutta et al. 2009; Mor et al. 2009; Alonso-Herrero et al. 
2011, 2013; Lira et al. 2013; Mori et al. 2014; Ruschel-Dutra 
et al. 2015; Ichikawa et al. 2015; Carcfa-Bernete et al. 2014). 
In this work we use an updated version of 


BAYESCLUMPY, a Bayesian tool specifically designed to fit 
the CLUMPY models to observed SEDs and mid-IR spec¬ 
troscopy (Asensio Ramos & Ramos Almeida 2009). Within 
BAYESCLUMPY there are two additional parameters not listed 
in Table [5] that can be fitted or fixed. The first one is the 
scaling factor called shift, which accounts for the vertical 
displacement needed to match the fluxes of a given model 
to an observed SED and/or spectrum. This is proportional 
to the AGN bolometric luminosity. The second one is the 
foreground extinction Ay due to the host galaxy which is 
different from that produced by the torus along the LOS. 
We use the Chiar & Tielens (2006) extinction law derived 
for the local ISM, as it provided a better fit to the IRS spec¬ 
trum (Section 3.2.1). Foreground extinction has been shown 
to be required to fit the nuclear IR emission of Seyfert nu¬ 
clei with deep silicate features (see e.g. Alonso-Herrero et 
al. 2011; Gonzalez-Martfn et al. 2013), as is the case for 
UGC 5101. We assume a uniform prior distribution for both 
of them, with the foreground extinction being in the range 
Ay (foreground) =0 — 40 mag. BAYESCLUMPY can also incor¬ 
porate the direct AGN emission to model the data, although 
in principle there is no clear detection of broad emission lines 
in UGC 5101 (see next section). 

Under the assumption that the torus models are valid 
simultaneously for photometric and spectroscopic data and 
that the noise in all observed data is uncorrelated, the 
Bayesian approach is insensitive to the fact that there are 
many more spectroscopic points than photometric ones. 
This is justified by the fact that the method samples the 
full posterior distributions and obtains marginalized pos¬ 
teriors for each parameter distribution (Asensio Ramos & 
Ramos Almeida 2009). The results of the fitting process are 
the marginal posterior distributions for the six free param¬ 
eters that describe the CLUMPY dusty torus model plus 
the vertical shift and foreground extinction. 

4.2 Results from the fit to the unresolved SED 
plus the starburst-subtracted spectroscopy 

We modelled the spectral energy distribution (SED) and 
spectroscopy of UGC 5101 using the near- and mid-IR un¬ 
resolved fluxes plus the starburst-subtracted GTC/CC spec¬ 
trum. We used as an upper limit the nuclear magnitude in 
the L-band from Imanishi et al. (2014). Both the observed 
unresolved SED and GTC/CC nuclear starburst-subtracted 
spectrum of UGC 5101 are plotted in Fig. [S] 

After fitting the unresolved emission and the starburst- 
subtracted nuclear spectrum, we obtained the marginal pos¬ 
terior distributions for the six free parameters that describe 
the CLUMPY dusty torus model (Table 5) plus the vertical 
shift and foreground extinction Ay. If the observational data 
contain sufficient information for the fit, then the resulting 
marginal posterior distributions are centred or show trends 
at certain values within the considered intervals. From these 
posterior distributions we also derived other torus parame¬ 
ters listed in Table 6. 

The posterior distributions of the fitted torus model pa¬ 
rameters show that they are well constrained (see Fig.Qand 
Table 5). Both the high value of the angular width of the 
torus (Ttorus ~ 60 degree and number of clouds along the 
equatorial direction No ^ 12 produce high apparent and ge¬ 
ometrical covering factors f{i) ~ 0.94 and /2 ~ 0.94, which 
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Table 5. Parameters of the CLUMPY torus models and for the best fit to nuclear emission of UGC 5101. 


Parameter 

Symbol 

Interval 

Median 

MAP 

Torus angular width [degrees] 

<^torus 

[15, 70] 

60lf 

66 

Torus radial thickness 

Y 

[5, 100] 


87 

Number of clouds along an equatorial ray 

No 

[1, 15] 

12tl 

13 

Index of the radial density profile 

<? 

[0, 3] 

o-stSI 

1.0 

Viewing angle [degrees] 

i 

[0, 90] 

55li« 

46 

Optical depth per single cloud 

TV 

[5, 150] 


32 


Notes.—Columns 1, 2 and 3 show the name of the CLUMPY torus model 
parameters, their symbol and the fitting range. Columns 4 and 5 show the median 
and maximum-a-posterior values obtained from the posterior distributions. 



Figure 6. Left: observed unresolved nuclear SED (red dots) and mid-IR spectroscopy after subtracting the starburst contribution from 
the CC nuclear spectrum (black dots). The orange dot is the nuclear 880/rm continuum flux from Wilson et al. (2008). We did not use 
this data point for the BayesClumpy fit (see text). The blue line and shaded region represent the best CLUMPY torus model and the 
range of models within the 68 per cent uncertainty in the best fitted parameters, whereas the solid orange line is the MAP model. Right: 
enlarged view of the best-fit CLUMPY torus models (as in the left panel) around the 9.7/rm and 18/rm silicate features. 


are defined as (see Alonso-Herrero et al. 2011 for details): 


fli) — ^-torusji) 

and, 


( 2 ) 


/2 = 1 — / Pesc{l3)cOs{(3)dl3. (3) 

Jo 

These high covering factors are consistent with the de¬ 
rived low escape probability for an AGN produced photon 
Pesc (see Table 6), which might explain the lack of strong 
AGN signatnres in the optical spectrum of UGC 5101. The 
htted angular size and number of clouds along the equato¬ 
rial direction of the UGC 5101 torus are consistent with the 
values found for Seyfert 2 nuclei, whereas the torus size of 
UGC 5101 35 pc) is larger than the typical sizes derived 

for Seyfert 1 and 2 nuclei using the same modeling tech¬ 
niques (Ramos Almeida et al. 2011, Alonso-Herrero et al. 
2011, Ichikawa et al. 2015). 

The high value of the foreground extinction, Ay = 


32jly mag, is necessary to reproduce the deep 9.7 fim sili¬ 
cate feature observed in the nuclear region. However, com¬ 
pared with the foreground extinctions derived from the 
Spitzer/IRS and CC spectral decompositions this value is 
nearly 3 times lower. This difference can be explained con¬ 
sidering that the CLUMPY models include the geometrical 
properties of the torus and, on this model, the foreground 
extinction is not the only contributor to the observed deep 
9.7 jj,m silicate feature. The best-fit models also reproduce 
well the 17.7 pm unresolved flux around the 18 pm silicate 
feature. All these values are also consistent with the param¬ 
eters, including the high covering factors, derived in other 
LIRGs and ULIRGs hosting deeply embedded AGN, e.g. 
Arp 299A (Alonso-Herrero et al. 2013), NGC 6240 (Mori et 
al. 2014), and Mrk 1066 (Ramos Almeida et al. 2014b). 

We translated the htted torus model parameters into 
two model spectra. The hrst one corresponds to the 
maximum-a-posterior (MAP) values that represent the best 
ht to the data. The second one is produced with the me¬ 
dian values of the probability distributions of the model 
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Figure 7. Marginal posterior distributions of the fitted CLUMPY torus model parameters to the nuclear emission of UGC 5101. The 
dotted lines represent the Icr intervals, whereas the solid black and orange lines are the MAP and median values of the distributions, 
respectively. 


parameters. We plot these best-fit CLUMPY models in 
Fig- [altogether with the both unresolved SED and starburst- 
subtracted GTC/CC spectrum. We also show in this figure 
the 880 pm continuum flux of 37 ± 9 mjy measured with a 
1 arcsec x 0.7 arcsec beam (Wilson et al. 2008), although we 
did not use it for the fit. Clearly, the torus emission is well 
below this flux, indicating that in the far-IR the nuclear 
emission is dominated by star formation activity, as already 
hinted by the spectral decomposition of the IRS spectrum 
(Section 3.2.1). Additionally, the 880pm flux may have a 
significant contribution from non-thermal or radio free-free 
emission (see discussion by Wilson et al. 2008). 

As can be seen from Fig. (6] the best-fit CLUMPY 
models reproduce well the shape of the silicate feature and 
the mid-IR unresolved emission. The 2.2 pm is not well fit¬ 
ted by the CLUMPY models. One possibility is that the 
near-IR unresolved emission is still contaminated by stellar 
emission from recent star formation in the nuclear region 
of the galaxy. However, in this case, the SED of the unre¬ 
solved emission should be bluer than the observed nuclear 
H — K = 1.9 color derived in this work, and at least 2 or¬ 
ders of magnitude more luminous than the AGN. Another 
explanation, more consistent with the nuclear H — K color, 
could be the presence of very hot dust clouds of graphite 
(T ~ 1000 K) (see for example, Mor & Trakhtenbrot 2011). 
In fact, Alonso-Herrero et al. (2011) found a slight excess at 
near-IR emission in the modelling of the SED+spectroscopy 
of some Seyfert Is and, Mor et al. (2009) also consider a hot 
dust component for modelling the SED-fspectroscopy of a 
sample of PG quasars. 

Thus, may be in UGC 5101, a combination of 
CLUMPY torus emission-l-cold foreground absorbing dust 


Table 6. Derived properties from the posterior distributions of 
the fitted CLUMPY torus parameters. 


Parameter 

Median 

Extinction-corrected Lboi [erg 

(1.9lo,?) X 10'*® 

Escape probability Pesc 

U.UUUO_Q QQQ5 

Geometrical covering factor /2 

u.y^_o.o7 

Apparent covering factor f{i) 

u.y^_o.i5 

Foreground extinction Ay [mag] 

32l^ 


Notes.— See equation 2 and 3 for the definition of the covering 
factors and escape probability. 

screen might not be the appropriate model for active nu¬ 
clei with very deep silicate features (e.g. Alonso-Herrero et 
al. 2011; Mori et al. 2014; Merlo et al. 2014) and spheri¬ 
cally symmetric smooth models may be more appropriate for 
deeply embedded objects (see Levenson et al. 2007; Alonso- 
Herrero et al. 2013). 

4.3 The AGN luminosity 

From the scaling of the fitted CLUMPY models to the 
data we derived an AGN bolometric luminosity of Lboi = 
(l.Oto r) 10"*^ erg s”*^. Using the hard X-ray vs. 12 pm 
correlation observed for local Seyfert galaxies (Gandhi et 
al. 2009; Levenson et al. 2009; Asmus et al. 2011) and 
the emission contribution by the starburst-subtracted nu¬ 
clear component at 12pm from the GTC/GC spectral de¬ 
composition after correcting for the derived foreground ex¬ 
tinction (see Table 4), we predict an intrinsic hard X-ray 
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luminosity L 2 _iokeV = (4.1 ± 0.7) x lO"'® erg s“^. Ap¬ 
plying a bolometric correction [10 — 20] x Z/ 2 -iokeV (Mar¬ 
coni et al. 2004) we found a bolometric luminosity between 
4.1 X 10"*^ ^ Lboi ^ 8.2 X 10“*^ erg s”*^. This range of AGN 
bolometric luminosities is lower than our estimate from the 
CLUMPY torus model fit to the unresolved nuclear IR emis¬ 
sion and the starburst-subtracted nuclear CC spectrum of 
UGC 5101. 

Gonzalez-Martm et al. (2009) on the other hand, es¬ 
timated an X-ray luminosity L 2 - lokev = 7 X 10"^^ erg s ^ 
after applying a Compton-thick correction to the observed 
2 — lOkeV emission. If we consider the above bolometric 
correction, then the AGN bolometric luminosity is con¬ 
sistent with the value derived from fitting the CLUMPY 
torus models to the unresolved IR emission and starburst- 
subtracted nuclear CC spectrum with that predict from the 
hard X-ray vs. 12 fim correlation. 

The inferred AGN bolometric luminosity from the 
CLUMPY torus modelling accounts approximately for ~ 
56% of the IR luminosity of the system, which is deter¬ 
mined using the fluxes in all four IRAS bands (12, 25, 60 
and 100 /rm) and the flux-luminosity correlation in Sanders 
& Mirabel (1996). Considering the uncertainties in the AGN 
bolometric luminosity (see Table 6) this value is consistent 
with other estimates in the literature. For instance, Lonsdale 
et al. (2003) using high angular resolution radio observations 
estimated a 30% AGN contribution in UGC 5101. Veilleux 
et al. (2009) derived a similar value (35%) using a number 
of spectral mid-IR diagnostics using the Spitzer/IRS spec¬ 
troscopy. 

In summary, despite the fact that the CLUMPY mod¬ 
els do not fit simultaneously all the near- and mid-IR unre¬ 
solved emission, and starburst-subtracted nuclear CC spec¬ 
trum of UGC 5101, the derived AGN bolometric luminos¬ 
ity is consistent with estimates using other methods. In our 
modelling the AGN bolometric luminosity is only a func¬ 
tion of the scaling factor of the models to the data, which 
is driven mostly by the fluxes near the peak of the torus 
emission. In the case of UGC 5101, the peak of the best-fit 
CLUMPY torus model is at Arest ~ 40 — 50 ^m (see Fig. 
left panel) and the Q-band data point (observed 17.7 /rm) is 
the closest to that peak. Therefore even for deeply embedded 
AGN, the modelling of the unresolved nuclear emission and 
the starburst-subtracted spectroscopy with the CLUMPY 
torus models can still be useful for deriving the AGN bolo¬ 
metric luminosity, as this quantity is driven by the scaling 
of the model and is less subject to the detailed modelling of 
the near and mid-IR SED and the silicate feature. 


5 SUMMARY AND CONCLUSIONS 

We have presented new high angular resolution (0.3 — 
0.4arcsec) 8.7/rm mid-IR imaging and 7.5 — 13 pm spec¬ 
troscopy of the nuclear region of the ULIRG UGC 5101 
using CC on the GTC. We also analyzed archival 
HST/mCMOS near-IR and Subaru/COMICS Q-band 
imaging and Spitzer/IKS spectroscopy. The presence of PAH 
emission in the nuclear (central ~ 393 pc) CC spectrum indi¬ 
cates that some of the dust in central region is being heated 
by the UV radiation of relatively young stars. Therefore, 
to study the physical and geometrical properties of the ob¬ 


scuring material surrounding the AGN of UGC 5101, it was 
necessary to disentangle the IR emission presumably due 
to dust heated by the AGN from that due to nuclear star 
formation and/or stellar emission. 

By using GALFIT we were able to separate the unre¬ 
solved nuclear emission from the extended host galaxy emis¬ 
sion using the near-IR (1.60 and 2.22 pm) and mid-IR (8.7 
and 17.7 pm) imaging data. We also calculated the unre¬ 
solved emission at 8.7 pm with other methods, namely PSF- 
subtraction and PSF-scaling and found that the resulting 
fluxes were similar to those derived with GALFIT. In addi¬ 
tion, we found that the Sersic parameters, which represent 
the extended emission, derived at near- and mid-IR wave¬ 
lengths are consistent with each other. 

For the mid-IR spectroscopy we used a spectral decom¬ 
position method to estimate the emission of the obscuring 
material around the AGN, which we represented as an ex¬ 
tinguished power-law continuum. We also allowed both com¬ 
ponents to have different foreground extinctions. From the 
spectral decomposition of the IRS data we found an index 
of the power-law (c 2 = 2.0 ± 0.2) consistent with the values 
observed for Seyfert 1 and Seyfert 2 nuclei (Ramos Almeida 
et al. 2011). Given the limited spectral range of the CC 
spectrum for the decomposition we fixed the index of the 
power-law and starburst template to those derived for the 
IRS data. Within the uncertainties the starburst-subtracted 
nuclear fluxes at 8.7 and 12 ^m CC continuum component 
are consistent with those derived from the spectral decom¬ 
position of the IRS spectrum. In addition, the starburst- 
subtracted nuclear flux at 8.7 pm is also consistent within 
the uncertainties with the value derived with the GALFIT 
modelling of the CC imaging data. We also estimated that in 
the central 393 pc of UGC 5101 the power-law component at 
8.7 pm and 12 pm (observed wavelengths) contributes 67% 
and 59% of the nuclear emission, respectively. 

Before fitting the unresolved SED of UGC 5101 with the 
CLUMPY torus models of Nenkovaet al. (2008a, b), we sub¬ 
tracted the fitted starburst component from the GTC/CC 
nuclear observed spectrum. We found that CLUMPY torus 
models with a high covering factor and a high foreground 
extinction (Av ~ 32 mag) fit well the silicate feature and un¬ 
resolved mid-IR emission. The high value of the foreground 
extinction is similar to results for other deeply embedded 
AGN in local LIRGs and ULIRGs (see e.g. Alonso-Herrero 
et al. 2013 and Mori et al. 2014). However, the observed 
unresolved near-IR emission of UGC 5101 is well above the 
fitted CLUMPY models. This could be due to contamina¬ 
tion by nuclear stellar emission even at the angular resolu¬ 
tion of O.lSarcsec of the NICMOS data or more likely due 
to the presence of hot dust clouds. We also showed that the 
observed 880 pm nuclear flux is much higher than the ex¬ 
trapolation to the far-IR of the fitted torus models. Again, 
this indicates that there is strong nuclear star formation, as 
also indicated by the presence of PAH emission in the CC 
spectrum, and this component also dominates the far-IR nu¬ 
clear emission of UGC 5101. 

Even though the CLUMPY torus models are not able 
to reproduce simultaneously the nuclear unresolved near and 
mid-IR emission, and the starburst-subtracted GTC/CC nu¬ 
clear spectrum of UGC 5101, the derived AGN bolometric 
luminosity (Lboi ~ 1.9 x 10“*® erg s”*^) is consistent with that 
estimated from hard X-rays (Gonzalez-Martm et al. 2009) 
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but higher than that from the 2 — lOkeV vs. 12 /rm correla¬ 
tion (Gandhi et al. 2009; Levenson et al. 2009; Asmus et al. 
2011 ). 


ACKNOWLEDGEMENT 

We thank the referee for a thorough report that has im¬ 
proved the paper significantly. This work has been partly 
supported by Mexican CONACyT under research grant 
CB-2011-01-167291. MM-P acknowledges support by the 
CONACyT PhD fellowship program. AA-H and AH-C ac¬ 
knowledge financial support from the Spanish Plan Na- 
cional de Astronomi'a y Astrofisica under grant AYA2012- 
31447, which is party funded by the FEDER programme, 
and the Universidad de Cantabria through the Augusto G. 
Linares programme. CRA is supported by a Marie Curie In- 
tra European Fellowship within the 7th European Commu¬ 
nity Framework Programme (PIEF-GA-2012-327934). IG- 
B ackowledges financial support from the Institute de As- 
trofisica de Canarias through Fundacion La Caixa and from 
the Spanish Ministry of Science and Innovation (MICINN) 
through project PN AYA2013-47742-c4-2-P (Estallidos). 

This work is based on observations made with the 
GTC, installed in the Spanish Observatorio del Roque de 
los Muchachos of the Institute de Astrofisica de Canarias, 
in the island of La Palma. It is also based partly on obser¬ 
vations obtained with the Spitzer Space Observatory, which 
is operated by JPL, Caltech, under NASA contract 1407. 
This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by JPL, Caltech, 
under contract with the National Aeronautics and Space 
Administration. The CASSIS is a product of the Infrared 
Science Center at Cornell University, supported by NASA 
and JPL. Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the data archive at 
the Space Telescope Science Institute. STScI is operated by 
the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS 5-26555. 


REFERENCES 

Allamandola, L. J., Hudgins, D. M., & Sandford, S. A., 1999, 
ApJL, 511, L115 

Alonso-Herrero, A., Quillen, A. C., Simpson, C., Efstathiou, A., 
& Ward, M. J., 2001, AJ, 121, 1369 

Alonso-Herrero, A., Ramos Almeida, C., Mason, R., et al., 2011, 
ApJ, 736, 82 

Alonso-Herrero, A., Pereira-Santaella, M., Rieke, G. H., & 
Rigopoulou, D., 2012, ApJ, 744, 2 

Alonso-Herrero, A., et al., 2013, ApJ, 779, L14 

Armus, L., Charmandaris, V., Spoon, H. W. W., et al., 2004, 
ApJS, 154, 178 

Armus, L., Charmandaris, V., Bernard-Salas, J., et al., 2007, 
ApJ, 656, 148 

Asensio Ramos, A., & Ramos Almeida, C., 2009, ApJ, 696, 2075 

Asmus, D., Gandhi, P., Smette, A., Honig, S. F., & Duschl, 
W. J., 2011, A&A, 536, AA36 

Asmus, D., Honig, S. F., Gandhi, P., Smette, A., & Duschl, 
W. J., 2014, MNRAS, 439, 1648 

Brandi, B. R., Bernard-Salas, J., Spoon, H. W. W., et al., 2006, 
ApJ, 653, 1129 


Buchanan, C. L., Kastner, J. H., Forrest, W. J., et ah, 2006, AJ, 
132, 1890 

Chiar, J. E., & Tielens, A. G. G. M., 2006, ApJ, 637, 774 
Cohen, J. G., 1999, AJ, 117, 2428 
Downes, D., & Solomon, P. M., 1998, ApJ, 507, 615 
Gandhi, P., Horst, H., Smette, A., et al., 2009, A&A, 502, 457 
Garcia-Bernete, I., Ramos Almeida, C., Acosta-Pulido, J. A., et 
ah, 2015, MNRAS, 449, 1309 

Genzel, R., Lutz, D., Sturm, E., et ah, 1998, ApJ, 498, 579 
Gerakines, P. A., Schutte, W. A., Greenberg, J. M., & van 
Dishoeck, E. F., 1995, A&A, 296, 810 
Gonzalez-Martin, O., Masegosa, J., Marquez, 1., & Guainazzi, 
M. 2009, ApJ, 704, 1570 

Gonzalez-Martin, O., Rodriguez-Espinosa, J. M., Diaz-Santos, 
T., et al. 2013, A&A, 553, AA35 
Haan, S., Surace, J. A., Armus, L., et ah, 2011, AJ, 141, 100 
Hernan-Caballero, A., Alonso-Herrero, A., Hatziminaoglou, E., 
et ah, 2015, ApJ, 803, 109 

Honig, S. F., Kishimoto, M., Gandhi, P., et ah, 2010, A&A, 515, 
AA23 

Houck, J. R., Charmandaris, V., Brandi, B. R., et ah, 2004, 
ApJS, 154, 211 

Ichikawa, K., Packham, C., Ramos Almeida, C., et ah, 2015, 
ApJ, 803, 57 

Imanishi, M., Dudley, C. C. & Maloney, P. R., 2001, ApJL, 558, 
L93 

Imanishi, M., Terashima, Y., Anabuki, N., & Nakagawa, T., 
2003, ApJL, 596, L167 

Imanishi, M., & Saito, Y., 2014, ApJ, 780, 106 
Kim, D.-C. 1995, Ph.D. thesis, Univ. Hawaii 
Kormendy, J., & Sanders, D. B., 1992, ApJL, 390, L53 
Lebouteiller, V., Barry, D. J., Spoon, H. W. W., et ah, 2011, 
ApJS, 196, 8 

Levenson, N. A., Sirocky, M. M., Hao, L., et ah, 2007, ApJL, 
654, L45 

Levenson, N. A., Radomski, J. T., Packham, C., et ah, 2009, 
ApJ, 703, 390 

Lira, P., Videla, L., Wu, Y., et ah, 2013, ApJ, 764, 159 
Lonsdale, C. J., Smith, H. E., & Lonsdale, C. J., 1995, ApJ, 438, 
632 

Lonsdale, C. J., Lonsdale, C. J., Smith, H. E., & Diamond, P. J., 
2003, ApJ, 592, 804 

Marconi, A., Risaliti, G., Gilli, R., et ah, 2004, MNRAS, 351, 
169 

Mason, R. E., Lopez-Rodriguez, E., Packham, C., et ah, 2012, 
AJ, 144, 11 

Medling, A. M., U, V., Guedes, J., et ah, 2014, ApJ, 784, 70 
Merlo, M. J., Perlman, E. S., Nikutta, R., et ah, 2014, ApJ, 788, 
6 

Mor, R., & Trakhtenbrot, B., 2011, ApJL, 737, LL36 
Mor, R., Netzer, H., & Elitzur, M., 2009, ApJ, 705, 298 
Mori, T. L, Imanishi, M., Alonso-Herrero, A., et ah, 2014, PASJ, 
66 , 93 

Mullaney, J. R., Alexander, D. M., Huynh, M., Goulding, A. D., 
& Frayer, D., 2010, MNRAS, 401, 995 
Mullaney, J. R., Alexander, D. M., Goulding, A. D., & Hickox, 
R. C., 2011, MNRAS, 414, 1082 
Nardini, E., Risaliti, G., Watabe, Y., Salvati, M., & Sani, E., 
2010, MNRAS, 405, 2505 

Nenkova, M., Ivezic, Z., & Elitzur, M., 2002, ApJL, 570, L9 
Nenkova, M., Sirocky, M. M., Ivezic, Z., & Elitzur, M., 2008a, 
ApJ, 685, 147 

Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, Z., & Elitzur, 
M., 2008b, ApJ, 685, 160 

Packham, C., Telesco, C. M., Hough, J. H., & Ftaclas, C., 2005, 
Revista Mexicana de Astronomia y Astrofisica Conference Se¬ 
ries, 24, 7 

Nikutta, R., Elitzur, M., & Lacy, M., 2009, ApJ, 707, 1550 



High angular resolution mid-IR emission in UGC 5101 15 


Peng, C. Y., 2002, AJ, 124, 294 

Pope, A., Bussmann, R. S., Dey, A., et al., 2008, ApJ, 689, 127 
Ptak, A., Heckman, T., Levenson, N. A., Weaver, K., & Strick¬ 
land, D., 2003, ApJ, 592, 782 
Puget, J. L., & Leger, A., 1989, ARA&;A, 27, 161 
Radomski, J. T., Pina, R. K., Packham, C., Telesco, C. M., & 
Tadhunter, C. N., 2002, ApJ, 566, 675 
Radomski, J. T., Pina, R. K., Packham, C., et ah, 2003, ApJ, 
587, 117 

Ramos Almeida, C., Levenson, N. A., Rodriguez Espinosa, J. M., 
et ah, 2009, ApJ, 702, 1127 

Ramos Almeida, C., Levenson, N. A., Alonso-Herrero, A., et al., 
2011, ApJ, 731, 92 

Ramos Almeida, C., Alonso-Herrero, A., Levenson, N. A., et ah, 
2014a, MNRAS, 439, 3847 

Ramos Almeida, C., Alonso-Herrero, A., Esquej, P. et ah, 2014b, 
MNRAS, 445, 1130 

Reach, W. T., Megeath, S. T., Cohen, M., et ah, 2005, PASP, 
117, 978 

Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., et ah, 2009, 
ApJ, 692, 556 

Roche, P. F., & Aitken, D. K., 1984, MNRAS, 208, 481 
Ruschel-Dutra, D., Pastoriza, M., Riffel, R., Sales, D. A., & 
Winge, C., 2014, MNRAS, 438, 3434 
Sanders, D. B., Soifer, B. T., Elias, J. H., et ah, 1988a, ApJ, 
325, 74 

Sanders, D. B., Soifer, B. T., Elias, J. H., Neugebauer, G., Sz 
Matthews, K., 1988b, ApJL, 328, L35 
Sanders, D. B., Mirabel, 1. F., 1996, ARA&A, 34, 749 
Sargsyan, L., Weedman, D., Lebouteiller, V., et ah, 2011, ApJ, 
730, 19 

Scoville, N. Z., Evans, A. S., Thompson, R., et ah, 2000, AJ, 
119, 991 

Soifer, B. T., Neugebauer, G., Matthews, K., et ah, 2000, AJ, 
119, 509 

Spoon, H. W. W., Keane, J. V., Tielens, A. G. G. M., et ah, 
2002, A&A, 385, 1022 

Stierwalt, S., Armus, L., Surace, J. A., et ah, 2013, ApJS, 206, 
1 

Sturm, E., Lutz, D., Tran, D., et ah, 2000, A&A, 358, 481 
Surace, J. A., Sanders, D. B., Sz Evans, A. S., 2000, ApJ, 529, 
170 

Telesco, C. M., Ciardi, D., French, J., et ah, 2003, SPIE, 4841, 
913 

Thompson, R. L, Storrie-Lombardi, L. J., Weymann, R. J., et 
ah, 1998, arXiv:astro-ph/9810285 
Valiante, E., Lutz, D., Sturm, E., Genzel, R., & Chapin, E. L., 
2009, ApJ, 701, 1814 

Veilleux, S., Kim, D.-C., Sanders, D. B., Mazzarella, J. M., & 
Soifer, B. T., 1995, ApJS, 98, 171 
Veilleux, S., Rupke, D. S., Kim, D.-C., et ah, 2009, ApJS, 182, 
628 

Wilson, C. D., Petitpas, G. R., lono, D., et ah, 2008, ApJS, 178, 
189 

Yuan, T.-T., Kewley, L. J., Sz Sanders, D. B., 2010, ApJ, 709, 
884 


